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space, the localization of the period doubling and halving allows one to predict the set 
of parameters (amplitudes and phase lag) where ionization occurs. 
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1. INTRODUCTION 

Among systems at the atomic level, one-electron systems are the most fundamental and 
have proven to be a source of numerous surprising discoveries |T] : The multiphoton 
ionization of hydrogen in a strong microwave field |2] is such a simple system with 
surprisingly complex dynamics. With the development of theory of chaos, its stochastic 
and diffusional nature was elucidated [3], and intense research activity in the last three 
decades has resulted in a rather complete understanding of this problem [H EJ E] . In 
the recent years, the attention has shifted from understanding the physical process 
to manipulating or controlling it [3, [8] . In this context control refers to tailoring the 
physical behavior of dynamical systems (which generically show chaotic dynamics) using 
"knobs" (i.e., suitable external parameters). Identifying such knobs and understanding 
the mechanism by which they affect the dynamics is the ultimate goal of such research. 

Bichromatic pulses 0, HO] are natural tools in atomic control research because 
they offer practical control parameters such as polarization, amplitudes and phases 
[nilI2l|T5linilIllISlinilIHllISll2nil2I]. It has been shown that it is possible to use the 
direction of transport in a ratchet by varying the phase lag in the bichromatic pulse [22J. 
In this manuscript we consider the ionization behavior of a one-dimensional hydrogen 
atom driven by a strong bichromatic linearly polarized electric field which is modeled 
by the following one- dimensional Hamiltonian (x > 0) in atomic units 



where the indices / and h refer to the low and high frequency modes with frequencies 
Ioj and hu, respectively. These two modes are frequency locked (and denoted h:l), i.e. 
I and h are integers. The experimental results on ionization probability obtained in 
Refs. [TTl [T9] for two cases, mode lockings 3:1 and 3:2, show two very distinct regimes 
as the phase lag (ft is varied : For 3:1, the ionization probability shows a plateau in 
phase located at a rather high value of ionization probability, and for 3:2, the ionization 
probability shows no such plateau but rather a small value of ionization probability. 
These results were confirmed by quantum calculations on the one-dimensional model (PQ). 

The purpose of this manuscript is to show that these experimental observations 
are qualitatively and quantitatively captured using a periodic orbit analysis, which 
reveals the classical bifurcations responsible for ionization. Our analysis also allows 
the prediction of ionization at other values of parameters without resorting to large 
numerical simulations. A part of the results presented in this manuscript were 
announced in a Fast Track Communication |23j . 

This paper is organized as follows: First, in Sec. [2j we summarize the residue 
analysis of periodic orbits [24]. In Sec. [3] and Sec. HI we consider two cases, the 3:1 mode 
locking and 3:2 one, respectively. These two cases show drastically different ionization 
behavior which we explain through bifurcations of selected periodic orbits. In these 
two cases, we described the bifurcations (if any) in the system. We also compute the 
bifurcation surface which is defined by the set of parameters where a change of linear 
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stability (associated with a period doubling or halving) has occurred in the chaotic sea. 
In Sec. [5] we discuss the generalities of residue curve behavior for other mode lockings. 

2. Residue analysis of periodic orbits 

First, we map Hamiltonian into action-angle variables (of the unperturbed system 
Fh = Ft = 0) such that the principal quantum number n is reflected by action J. The 
relationship between J and n reads 

J = (t u) 1/3 n, 

where t = 2.4188843243 x i(T 17 s (atomic time unit) and u = 2nf = 12vr GHz in 
both cases. The action- angle variables [25] are denoted (J, 9) and obtained through the 
canonical change of coordinates x = 2 J 2 sin 2 ip, p = J -1 cot ip with 6 = 2ip — sin 2(p. 

We assume to = I without loss of generality (after a rescaling H' = cu~ 2 ^ 3 H, 
t' = tot, x' = to 2 ^ 3 x, p' = uj~ l l 3 p and consequently 4>' = 4> and F' = u~ 4 ^ 3 F). The 
Hamiltonian ([1]) becomes [I] 

H= -^ + 2J 2 [F h sm(ht) 

+ F t sin(Zt + <f))} (^a /2 + a k cos kO^j , (2) 

where a n = [J n (n) — J n _i(n)]/n and J n 's are Bessel functions of the first kind. Note 
that for a given mode locking h:l, there are three variables (J, 9, t) and three parameters 
(Fh, Ft, 4>). We denote this Hamiltonian H(J, 6, t; Fh, Ft, 4>). In general, for a given set of 
parameter values, the phase space as depicted on a Poincare section (a stroboscopic plot 
with period 2tt, see Fig. CD), is composed of a mixture of regular structures surrounded 
by a chaotic sea. More precisely, the lower part of phase space (J small) is composed of 
rotational invariant tori (which persist from the integrable case F% = Ft = as asserted 
by KAM theorem). The upper part is formed by regular islands surrounded by an 
unbounded chaotic sea. The ionizing trajectories are the ones in the chaotic sea which 
are unbounded (J, or equivalently n, becomes progressively large). At the center of 
the regular islands, there are elliptic (stable) periodic orbits which organize the regular 
motion around them. These periodic orbits partly result from the break-up of resonant 
tori into a pair of elliptic/hyperbolic orbits (according to Birkhoff's theorem). Other 
periodic orbits result from the bifurcation of these orbits. 

The general idea which is applied here is to follow a finite set of periodic orbits (both 
elliptic or hyperbolic) which have been identified as important. The criteria of choice 
combine several factors: the size of the island, the period and the location. For each 
periodic orbit of this set, we compute its location (J, 9) and its linear stability property 
as given by the residue (which is to be defined below). As the three parameters are 
varied, we follow the locations and residues of these orbits instead of computing the 
Poincare section for each value of parameters. This allows us to have a clear idea of 
what is going on in phase space and to predict ionization thresholds. 
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In what follows, we perform two kinds of computations : First, we compute the 
residue curves which are obtained as functions of the relative phase <ft f° r fixed values 
of the amplitudes (Fh,Fi). Second, we compute bifurcation surfaces which are defined 
as the set of parameters (F h , F h 0) where a change of linear stability has occurred. 

In order to start monitoring the stability of a family of periodic orbits, we first 
consider a specific periodic orbit, denoted O(0), of Hamiltonian (J2J) for <fi = which 
is our reference case. Numerically it is determined using a modified Newton-Raphson 
multi-shooting algorithm as described in Ref. [26]. The initial condition for launching 
the iterative algorithm can be taken from a Poincare section, for instance. As is 
continuously varied, the orbit O(0) deforms continuously into 0(<p), whose period is 
denoted T(<p). In addition to its location, we also monitor its linear stability properties 
given by the integration of the reduced tangent flow 

¥± = w*H{j,e,ti<f>)J}, 

where J = ^ ^ J ^ and V 2 H is the two-dimensional Hessian matrix (composed of 

second derivatives of H with respect to its canonical variables J and 9). The initial 
condition is J® = I2 (the two-dimensional identity matrix). The two eigenvalues of 
the monodromy matrix which make a pair (A(0), 1/A(0)) determine the stability 

properties (the flow is volume preserving, thus the determinant of is equal to 

1). The periodic orbit is elliptic if the spectrum is (e 1 "^, e - *^^) (stable, except for 
some particular cases), or hyperbolic if the spectrum is (A(0), 1/A(0)) with A(0) G M* 
(unstable). The stability properties can be deduced in a concise form using Greene's 
residue R [271 EE] 

9 - tr 7 T{4i) 

m = 

If R(4>) G]0, 1[, the periodic orbit is elliptic; if R((p) < or R((j)) > 1 it is hyperbolic; 
and if R(<fi) = and R(4>) = 1, it is parabolic. Generically, periodic orbits and their 
linear stabilities are robust against small changes of parameters, except at specific 
values where bifurcations occur [29]. These rare events affect the dynamical behavior 
drastically. They can be associated with an enhancement as well as a reduction of 
stability depending on the type of bifurcations, as shown in Refs. [211 EH]. We identify 
the bifurcations (if any) of a set of short periodic orbits, i.e., the type and the value of 
the parameter <fr c where they bifurcate. This provides a way to foretell if a relatively high 
ionization rate should be expected or not. The importance of considering two associated 
Birkhoff periodic orbits (i.e. periodic orbits with the same action but different angles in 
the integrable case, one elliptic and one hyperbolic), was emphasized in Ref. [23]. The 
main reason is that it allows one to discard some specific bifurcations (like collisions or 
exchanges of stability) and draw appropriate conclusions concerning the enhancement 
or reduction of stability. 
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Figure 1. Poincare section of Hamiltonian @ for Case (I) at <fi = 0. Full big circles 
(or big crosses, respectively ) indicate the two elliptic (resp. hyperbolic) periodic orbits 
with period 2ir we consider. Full small circles ( or small crosses, respectively) indicate 
the two elliptic (resp. hyperbolic) periodic orbits with period An we consider. The 
horizontal line corresponds to the principal quantum number n = 51. 

In the following sections, we analyze the short-time dynamics through residues of 
selected periodic orbits. We compute residue curves h- > R(4>, F h , F{) for fixed values of 
the amplitudes Fh and Fi, and bifurcation surfaces associated with a given periodic orbit 
defined as the set of parameters such that R(<p, F h , F{) = 1. We correlate these results 
with ionization probabilities obtained experimentally and also by quantum simulations 
through an empirical formula. 

3. 3:1 mode locking 

For this mode locking, we consider Fh = 24 Vcm -1 and Fi = 53.4 Vcm -1 as in Ref. [19|. 
These values correspond to the dimensionless values Fh = 24uo(io^) -4 ' 3 / do = 0.5278 and 
F l = 53.4^0 (t w)- 4/3 /d = 1-1743, where v = 0.036749326 and d = 1.88972613 x 10 s 
are both conversion factors, t = 2.4188843243 x l(T 17 s (atomic time unit). These 
values are referred as Case (J) in what follows. 

3.1. Poincare section 

Figure [T] shows a Poincare section of Hamiltonian (j2j) for Case (J) at = 0. The 
phase space is divided into two main parts : a lower regular containing many invariant 
tori, and a upper chaotic sea where trajectories escape rapidly to unbounded actions J 
(ionized trajectories). We notice two main islands in the chaotic sea. At the centers 
of these islands sit elliptic periodic orbits with period 2n (indicated by full circles). In 
addition, there is also a period 2 island in between these two main islands (associated 
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Figure 2. Poincare section of Hamiltonian @ for Case (I) at <fr — tt/3. Full circle 
(respectively cross) indicates the elliptic (resp. hyperbolic) periodic orbit with period 
27r we consider. The horizontal line corresponds to the principal quantum number 
n = 51. 



with an elliptic and hyperbolic periodic orbits with period Att). In the region of phase 
space around Jj ~ 0.49 (indicated by a straight line) where the initial states n = 51 
are prepared [19], any regular structures like the main islands and smaller ones are 
associated with trappings and hence reduce the ionization rate. 

Figure [2] shows a Poincare section of Hamiltonian (T5]) for Case (/) at = n/3. 
Apparently the upper elliptic periodic orbit of period 2ir in Fig. [1] has disappeared in 
chaotic sea. This case is associated with a higher ionization probability. Therefore, the 
loss of stability of the main periodic orbit near Jj plays a major role in the ionization 
process. 

The residue method will monitor the location and stability of elliptic periodic orbits 
as well as their associated hyperbolic periodic orbits [21] as a function of the parameters 
of the system. 

For Case (J), Figure [3] shows the positions of the upper elliptic and hyperbolic 
periodic orbits of period 2ir on the Poincare section as functions of 0. We notice that 
the action J is changing weakly as is varied. In contrast, its angle 9 is very sensitive 
to this parameter. 

3.2. Residue curve 

We follow the residue for each of the elliptic and hyperbolic periodic orbits mentioned 
above as the parameter <p varies for Case (J). Figure H] shows the residue curves of these 
orbits. 

In Fig. HI we monitor the upper elliptic periodic orbit (period 2tt) of Fig. [1] from 
= 0. This periodic orbit remains elliptic (R(4>) G]0, 1[) until C ~ 0.49 where a 
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Figure 3. The positions of the upper elliptic and hyperbolic periodic orbits with 
period 2ir on the Poincare section as functions of <f> for Case (7). The solid curves 
and the dashed ones correspond to the upper elliptic and hyperbolic periodic orbits of 
Fig. [1] respectively. 
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0.5 1 1.5 2 

Figure 4. Residue curves for the four periodic orbits with period 2n (solid curves) 
and the two periodic orbits with period 4-7T (dashed bold curves), indicated by crosses 
and circles on Fig. [T] for Case (I). The solid bold curves are for the upper set 
of elliptic/hyperbolic orbits of period 2~k. Small arrows indicate where bifurcations 
happen. The dotted bold curve between two arrows is associated with the residues of 
the elliptic periodic orbit with period 47r born of the period doubling bifurcation. 
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Figure 5. Bifurcation diagram for Case (J) showing the bifurcations indicated by 
arrows in Fig. 3J 



bifurcation occurs. At this critical point, the orbit turns parabolic. Increasing 
further, it turns and remains hyperbolic (R(4>) > 1) until 2ti/3 — C where another 
bifurcation appears, making the orbit elliptic again. This bifurcation process is of the 
period doubling kind at C and a period halving at 27r/3 — C (by symmetry). 

The bifurcation diagram for Case (J) is shown in Fig. [51 The computation of 
log | A-j- (<^>) | (where A±(0) are the two eigenvalues of the monodromy matrix 
associated with the upper elliptic periodic orbit) gives log|A±(0)| = before the 
bifurcation and log |A±(0)| oc ±a/0 — C right after the bifurcation. 

Figure [6] shows a projection of the upper elliptic periodic orbit undergoing the 
period doubling bifurcation in x — p representation where we note the doubling of the 
number of branches. 

Note that while the upper elliptic periodic orbit undergoes a bifurcation as is 
varied, the other three periodic orbits with period 2ir retain the stability properties 
they had at = 0. 

In the parameter range G]0 C , 2n/3 — C [, the upper part of phase space does 
exhibit more chaos as it can be shown on the Poincare section for = tt/3 (see Fig. [2]). 

Since the initial atomic beam is taken in the region with principal quantum number 
n — 51 which corresponds to action Jj w 0.495, the ionization rate is expected to 
be higher in the regime where there are no big islands in the chaotic sea, i.e. for 
G]0 c ,27r/3 — C [. On a finer scale, one has to take into account the smaller regular 
structures that are present in the chaotic sea, like for instance the period-2 island on 
the Poincare section (period Air). A similar period doubling and halving occur at 
Cj i ~ 0.199 and 2n/3 — Cj i ~ 1.895 respectively, as can be seen on Fig. 0] (upper 
bold dashed curve). We notice that the residue curve for the period Att orbit is higher 
than the upper solid curve (for the period 2tc orbit) and the phase region between period 
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Figure 6. Upper elliptic periodic orbit (period 2tt) of Case (I) undergoing the periodic 
doubling bifurcation in x — p representation. The insets shows the doubled number of 
branches indicating periodic doubling bifurcation. 



doubling (0 c ,i) and period halving (27r/3 — Cj i) is much wider than the one based on 
the period 2tt. This observation indicates that for longer pulse duration for which longer 
periodic orbits have to be taken into account, the period 47r orbit obtained by a repetition 
of the period 2tt orbit is more important (as a limiting factor) for ionization than the 
orbit with primary period 4n. As a by-product, the bifurcation at Cj i does not affect 
significantly the leading ionization behavior based on period 2tt. Another important 
periodic orbit with period 4ir is the one born of the period doubling bifurcation of the 
upper periodic orbit at <p c . This orbit is close to the newly hyperbolic periodic orbit 
and is at first elliptic, as seen on Fig. HI It is expected that the associated elliptic 
islands slow down ionization. This orbit experiences a period doubling bifurcation at 
0c,2 ~ 0.643 and halving bifurcation at 2tt/3 — <p c ^ ~ 1.451. In summary, the shortest 
periodic orbits (period 2n) always play the most dominant role on ionizations rates, 
although with larger pulse durations longer periodic orbits might be taken into account 
for a more detailed analysis. 

For values of <fi around 7r/3, a plateau is expected in the ionization probability 
versus <fi. The reason is that a strongly hyperbolic orbit only influences the ionization 
time and not the value of the ionization probability. Of course, this is true provided 
that the duration of the maximum pulse envelope is large enough. In the experiment, 
this is approximately 15 times the period of the shortest periodic orbits considered. 
Roughly speaking, this means that in the chaotic region all the orbits ionize (i.e., escape 
to a value of the action J ion > 1.26) regardless of the hyperbolicity degree. In Ref. [19J, 
experimental results as well as one-dimensional quantum calculations show this plateau. 
From quantum calculations, <p c ~ 0.5 was obtained in Ref. [19] which is in very good 
agreement with the parameter value <p c m 0.49 at which the bifurcation of the upper 
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Figure 7. Normalized ionization probability vs <j> based on Eq. ([3]) for Case (I) with 
A = —2.52 and B — 0.65. Circles represent the data obtained by one-dimensional 
quantum calculations, taken from Ref. [19j . Only periodic orbits with period 2n are 
considered. 



elliptic periodic orbit (period 2n) occurs. If the duration of the experiment or simulations 
is longer (two or three times longer), then higher-order regular structures (like the regular 
island of period Air in the Poincare section and the elliptic periodic orbit with period 
47r born of bifurcations) will play a role and we expect a similar, but smaller, plateau 
for g]0 Ci 2,2tt/3 - C>2 [. 

The periodic orbit analysis above elucidates whether or not there is a significant 
ionization probability for specific parameter values, and also where plateaus are expected 
to occur. This qualitative agreement highlights the important role played by these 
orbits. Furthermore, we can obtain quantitative agreement concerning the shape of the 
ionization curve versus phase lag by using the residue curves. Here we devise an 
empirical formula for relative ionization probability in the following way : First, the 
values of <ft giving the highest ionization would be the ones associated with the highest 
variations of the residues (in absolute value) with respect to the minimum ionization. 
Second, if the periodic orbit is too far (in action) from the considered action J; then it 
will not influence the dynamics so there should be a penalizing term depending on its 
position with respect to the chosen rescaled action. The relative ionization probability 
formula reads : 

p fx\ a i d ST- exp|fl m (0) -R m (4>o)\ ,^ 
P ion (0) = A + B ^ ===== , (3) 

m=1 exp|J m (0) - Ji\ 

where the sum is taken over the M different periodic orbits considered and J m (0) = 
J J (6)d6/27i is the action of the periodic orbit m. The parameters A and B in Eq. ([3]) 
are merely a translation and a dilatation of the curve in order to match the mean value 
and the amplitude of variations of P lon obtained in Ref. [TH]. This formula takes into 
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whereas the dashed one is for cf> = 1.95. 

account the value of the residues at 0o, where the minimum ionization takes place for 
each case. The aim is to set up a baseline for each of the periodic orbits (which is taken 
here at the value of the parameter where the ionization is minimal). Specifically for 
Case (/) 0o = [19]. In general, Eq. ([3]) can exhibit values which are greater than 1, 
which are not relevant. In order to remedy to this problem, we truncate Pj on at the value 
where a bifurcation occurs in accordance with the previous discussion on the relevance 
of the degree of hyperbolicity. Therefore in the range where R n {4>) is larger than one, 
Pion is constant (taken as the value of the residue at C where the bifurcation occurs). 

Figure [7J depicts P ion given by Eq. ([3]) versus parameter as well as the data 
taken from Ref. [19] for Case (/). Only the shortest periodic orbits with period 2-7T 
are considered as relevant. We notice that the empirical formula reproduces accurately 
the results obtained from quantum calculations. If we take into account the period 4-7T 
orbits for longer pulse duration, then a second plateau appears for &)<fic,2, 2tc/3 — 4> c ^[ 
together with an increase at Cj i and a decrease at 2n/3 — (f> c> \. 

3.3. Bifurcation surface 

The residue method is carried out to also predict the behavior of the system as all three 
parameters (the two amplitudes F h and Fi and the phase lag 0) are varied. In Fig. [HJ we 
represent the set of parameters where the upper elliptic periodic orbit (with period 2n) 
of Fig. [I] is in fact parabolic (i.e., the set of parameters where the system undergoes a 
major bifurcation). The equation of this surface in parameter space is P(0, Fh, Fi) = 1. 
The boundaries of the plateaus in parameter of Fig. [7J obtained by fixing the two 
values for Fh and Fi are on this surface. When Fh approaches zero, this surface is 
less dependent on parameter 0. Table [U reports some values based on our analysis 
which are in good agreement with experimental results from Ref. [T7j. Of course, the 
surface of Fig. [8] could also have been obtained from tedious numerical simulations of a 
large number of classical trajectories for each value of the parameters (0, P^P;). This 
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Table 1. Ionization thresholds obtained for Fh — 6 Vcm , experimentally in Ref. p"7] 
and by the residue method (see Fig. 2]). The If case corresponds to Fh = 0. 
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Figure 9. Poincare section of Hamiltonian (J2]) for Case (77) at cf> = 0. Full circles 
(respectively crosses) indicate the two elliptic (resp. hyperbolic) periodic orbits with 
period 2ir we consider. Small crosses indicate the two hyperbolic periodic orbits with 
period Air we consider. The horizontal line corresponds to the principal quantum 
number n = 51. 



integration needs to be performed for a sufficiently long time in order to decide if a 
given trajectory leads to ionization or not. In contrast, only one orbit for a short time 
(typically the period of the field) is needed for the residue analysis. Furthermore, using 
residues, this surface can be constructed locally without any need to consider all possible 
values of the parameters. 

4. 3:2 mode locking 

In what follows, we study in detail one set of values for the amplitudes of the fields 
Fh = 25 Vcm -1 and Fi = 33.5 Vcm -1 which correspond to Fh = 0.5498 and F\ = 0.7367 
in dimensionless units. This case will be referred as Case (//). 

4-1. Poincare section 

Figure [9] shows a Poincare section of Hamiltonian (j2J) for Case (II) at <fi = 0. Similar 
to Case (/) , there are two primary islands in the chaotic sea where two elliptic periodic 
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Figure 10. Residue curves for the four periodic orbits with period 2n indicated 
by crosses and circles on Fig. |H] and the two periodic orbits indicated by small 
crosses with period 4ir for Case (17). The bold solid curves are for the upper set 
of elliptic/hyperbolic orbits with period 2n. The thin solid curves are for the lower 
set of elliptic/hyperbolic orbits with period 2ir. The dashed curves are for the two 
initially hyperbolic periodic orbits with period 4tt. The inset shows irregular behavior 
of residue at 4> ~ 7r/6. The same behavior also occurs at 4> — 5tt/6 for that residue 
curve. 



orbits with period 2w sit at the centers, and the two associated hyperbolic orbits are 
in the chaotic sea. In addition, there are also two associated hyperbolic periodic orbits 
with period An as indicated by small crosses. One should notice that the (rescaled) 
principal quantum number considered in Ref. [1_9] lies in between these two islands (see 
the continuous horizontal line). We monitor the stability of this set of periodic orbits 
as we have performed for Case (J). 

4-2. Residue curve 

Figure [10] shows the four residue curves of period 2ti orbits and the two residue curves of 
period 4tt orbits. For short periodic orbits (period 2ir), the elliptic periodic orbits remain 
elliptic and the hyperbolic ones remain hyperbolic for all values of <fi. No bifurcation 
occurs except in a small range of phase (see inset of Fig. [TOl) where no significant 
stability change is observed. Consequently, the ionization probability is expected to 
be approximately independent of (f> and to be lower than Case (I) since for these values 
of amplitudes, the chaotic region is smaller. This is consistent with the experimental and 
quantum calculations of Ref. |19j . The experimental results show a nearly flat curve for 
the ionization probability versus 0, whereas the quantum calculations show significant 
variations for this probability but no sharp increase and decrease as in Case (/). 

When the duration of pulse is longer we should take into account the effect of longer 
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Figure 11. Normalized ionization probability vs <fi based on Eq. for Case (II) with 
A = —0.485 and B — 0.17. Circles represent the data obtained by one-dimensional 
quantum calculations, taken from Ref. [19j . Only periodic orbits with period 2ir are 
considered. 

periodic orbits whose residues are shown by bold dashed curves in Fig. [10] with period 
4-7T in this particular case. Apparently the upper period Air orbit experiences period 
halving bifurcations at ~ 0.07 and ~ 2.164; and period doubling bifurcations at 
w 0.977 and « 3.072. 

Figure [TT1 depicts P ion given by Eq. (j3J) versus parameter as well as the data taken 
from Ref. [19] for Case (II). Only periodic orbits with period 2ir are considered. In 
Eq. ([3]) we again take 0o — since the ionization is minimal at = for Case (II) 
|19j . We notice that it captures some essential features of the ionization curve, like 
the two unequal-sized peaks and the specific shape of both peaks (one more peaked, 
the other one, more round). This feature results from the asymmetry property of 
bichromatic microwave amplitude [30J. For this case, the bichromatic microwave field at 
= 7r/6±A0 has exactly the same amplitude but opposite direction as at = 7r/2±A0, 
where < A0 < ir/6. Replacing "+" in front of the bichromatic field by "-" in Eq. (TjQ) 
leads to a Hamiltonian for the other direction. Because of this symmetry property of 
bichromatic field, the ionization probability along the other direction can be obtained 
from Fig. [11] by a horizontal translation of tt/3. The total ionization rate should be the 
sum of the rates for both directions, which shows equally high peaks for total ionization 
rate in 0. For longer pulse duration, we should consider longer periodic orbits like those 
with period Air. Again we should take the baseline (minimal ionization point) at = 0. 
Since in Fig. [10] the residue value of elliptic periodic orbit with period An (upper dashed 
curve) at minimal ionization point (0 = 0) is already greater than 1, the ionization 
rate should be high and roughly identical on the whole space. This indicates that 
when pulse duration gets longer, the ionization rate for any goes higher and suggests 




Periodic orbit bifurcations as an ionization mechanism 



15 



ionization rate for any <fi will reach a roughly constant upper limit if the time duration is 
sufficiently long. This property is obviously true because when the system is exposed to 
external microwave field for a longer time, more electrons will gain sufficient energy to 
escape. However, as mentioned in Sec. [3] the shortest periodic orbits (period 2ti) always 
play the most dominant role on ionizations rates, while the longer periodic orbits might 
be taken into account heuristically for a finer analysis when the pulse duration gets a 
lot longer. Finally, our analysis also indicates the asymmetric ionization property does 
not exist for 3:1 mode locking case, in agreement with Ref. [30J. 

5. Generalization to h:l mode locking 

In order to generalize this approach, we also investigate h:l— 2:1, denoted Case (III), 
and 5:1, denoted Case (IV) mode locking. 

5.1. Residue curves 

For Case (III), we take = 24 Vcm -1 , Fi = 53.4 Vein -1 and the high frequency of 
12 GHz. Figure [12] shows Poincare section at = 0. It indicates the relevant periodic 
orbits to consider in the residue analysis. Figure [13] shows the residue curves based on 
the two periodic orbits with period 2n in the chaotic sea of Fig. [12l No bifurcations 
take place for these orbits. Figure [TH depicts relative ionization rate with respect to <j) 
based on Eq. ([3]). For this case we take 0o — in Eq. ([3]) since the minimal ionization 
point is at <ft = according to the maximum field rule (see next section). Based on our 
approach, the asymmetric ionization property along a single direction (unequal-sized 
peaks) does appear for the 2 : 1 mode locking, in agreement with Ref. [30] . 

For Case (IV), we take Fh = 24 Vcm -1 , Fi = 53.4 Vcm -1 and the high frequency 
is 30 GHz. Figure [T5l shows a Poincare section at <fi = 0, and Fig. [TBI shows the residue 
curves based on the six periodic orbits of period 2tt in the chaotic sea of Fig. [15] We 
notice that the residue based on the upper elliptic periodic orbit has similar behavior to 
the one of the lower elliptic periodic orbit for Case (II) as shown in the inset of Fig. [101 
The two residue curves (two thin lines) in Fig. [16] from the two lower periodic orbits in 
Fig. [15] are almost constant and therefore have little influence on ionization probability. 
Figure [T7] depicts relative ionization rate with respect to <fi based on Eq. ([3]) according 
to Fig. [16] (The lowest two periodic orbits are not considered since their corresponding 
residue curves are almost constant). For this case we take 0o = tt/5 in Eq. (J3j) since 
the minimal ionization point is at <fi = it/ 5 according to the maximum field rule. The 
asymmetry ionization property does not exist for the 5:1 mode locking. 

5.2. Comparison with the maximum field rule 

A traditional way of explaining ionization probability behavior is the maximum field 
rule (peak field amplitude rule) formulated in Ref. [J71 [3TJ . The relative ionization 
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Figure 12. Poincare section of Hamiltonian for Case at <p = 0. Full circle 

(respectively cross) indicates the elliptic (resp. hyperbolic) periodic orbit with period 
27r we consider. The horizontal line corresponds to the principal quantum number 
n = 51. 
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Figure 13. Residue curves for the two periodic orbits with period 2ir indicated by 
cross and circle on Fig. [12] for Case (III). 

probability vs <fi based on the maximum field rule can be written in the form 

ikm(</>) = C + Drasx^ [QM \F h wn.{ht) + F L sm(lt + <j>)\. (4) 

The parameters C and D in Eq. (j3J) are merely a translation and a dilatation of the 
curve in order to make the curves from the maximum field rule comparable to those 
based on bifurcation analysis. Figure [18] shows relative ionization probability based on 
Eq. (jlj) for four different cases. Compared with our periodic orbit bifurcation analysis in 
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Figure 14. Relative ionization probability vs cf> at principal quantum number n = 51 
based on Eq. © for Case (J//) with A = and B = 0.1. The absolute ionization rate 
may differ according to the pulse duration. 
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Figure 15. Poincare section of Hamiltonian @ for Case (IV) at 4> = 0. Full circles 
(respectively crosses) indicate the elliptic (resp. hyperbolic) periodic orbits with period 
2n we consider. The horizontal line corresponds to the principal quantum number 
n = 47. 
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Figure 16. Residue curves of the six periodic orbits with period 2tt indicated by 
crosses and circles on Fig. [T5] for Case (IV). 
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Figure 17. Relative ionization probability vs </> at principal quantum number n = 47 
based on Eq. © according to Fig.[16]for Case (IV) with A = -0.2 and B = 0.1. The 
absolute ionization rate may differ according to the pulse duration. 

finite pulse duration, apparently the maximum field rule does not result in quantitative 
agreement with the quantum simulation results of Ref. [19] for either Case (/) or (//). 
In Case , the maximum field rule does not produce the two unequal-sized peaks as 
our bifurcation analysis does. Only in Case (IV) do the results based on the maximum 
field rule agree with our bifurcation analysis because there is neither bifurcation nor 
asymmetric ionization property for this case. Insets of Fig. [T8lshow one cycle (t e [0, 2ir]) 
bichromatic field e(t) = Fhsin(ht) + -F;sin(/t + <fi) with typical parameters as used for 
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our previous analyses in dimensionless units for different values. Arrows indicate the 
values for which the bichromatic fields are drawn. Generally for Case (I) the positive 
component of the bichromatic field is the same as the negative component for any 
in each full cycle (directional symmetry), whereas this directional symmetry does not 
appear for Case (III) [30] except for some specific values like = 0,1, n. Similarly, for 
Case (II) this directional symmetry generally does not appear except for some specific 
values like — 0, |, y, 7r, -y, whereas for Case (IV), this directional symmetry does 
appear for any values. Although the maximum field rule can be used to determine 
qualitative features like the minimal ionization point with respect to 0, the quantitative 
agreement is not as satisfactory as the one given by a method which relies on analyzing 
the chaotic dynamics like the one used in this article based on periodic orbit bifurcation 
analysis to depict relative ionization. 

Acknowledgements 

This research was supported by the US National Science Foundation. C.C. acknowledges 
support from Euratom-CEA (contract EUR 344-88-1 FUA F). We thank Xavier Leoncini 
and Luca Perotti for useful discussions. 

References 

[1] R. Bliimel and W.P. Rcinhardt, Chaos in Atomic Physics (Cambridge University Press, Cambridge, 
UK, 1997). 

[2] J.E. Bayfield and P.M. Koch, Phys. Rev. Lett. 33, 258 (1974). 

[3] B.I. Meerson, E.A. Oks and P.V. Sasorov, Pis'ma Zh. Eksp. Teor. Fiz. 29, 79 (1979) [JETP Lett. 
29, 72 (1979)]. 

[4] G. Casati, B.V. Chirikov, I. Guarneri and D.L. Shepelyansky, Phys. Rep. 154, 77 (1987). 

[5] R.V. Jensen, S.M. Susskind and M.M. Sanders, Phys. Rep. 201, 1 (1991). 

[6] P.M. Koch and K.A.H. van Leeuwcn, Phys. Rep. 255, 289 (1995). 

[7] H. Rabitz, R. de Vivie-Riedle, M. Motzkus and K. Kompa, Science 288, 824 (2000). 

[8] M. Shapiro and P. Brumer, Phys. Rep. 425, 195 (2006). 

[9] L. Ko, M.W. Noel and T.F. Gallagher, J. Phys. B 32, 3469 (1999). 
[10] F. Ehlotzky, Phys. Rep. 345, 175 (2001). 
[11] J.E. Howard, Phys. Lett. A 156, 286 (1991). 

[12] A. Haffmans, R. Bliimel, P.M. Koch and L. Sirko, Phys. Rev. Lett. 73, 248 (1994). 

[13] A. Buchleitner, D. Delandc and J.-C. Gay, J. Opt. Soc. Am. B 12, 505 (1995). 

[14] D. Petrosyan and P. Lambropoulos, Phys. Rev. Lett. 85, 1843 (2000). 

[15] M. Ivanov, P.B. Corkum, T. Zuo and A.D. Bandrauk, Phys. Rev. Lett. 74, 2933 (1995). 

[16] L. Sirko, S.A. Zelazny, and P.M. Koch, Phys. Rev. Lett. 87, 043002 (2001). 

[17] L. Sirko and P.M. Koch, Phys. Rev. Lett. 89, 274101 (2002). 

[18] V. Batista and P. Brumer, Phys. Rev. Lett. 89, 143201 (2002). 

[19] P.M. Koch, S.A. Zelazny and L. Sirko, J. Phys. B 36, 4755 (2003). 

[20] C. Rangan, A.M. Bloch, C. Monroe and P.H. Bucksbaum, Phys. Rev. Lett. 92, 113004 (2004). 
[21] V. Constantoudis and C.A. Nicolaides, J. Chem. Phys. 122, 084118 (2005). 

[22] G. G. Carlo, G. Benenti, G. Casati and D. L. Shepelyansky, Phys. Rev. Lett. 94, 164101 (2005). 

[23] S. Huang, C. Chandre and T. Uzer, J. Phys. B 40, F181 (2007) 

[24] R. Bachelard, C. Chandre and X. Leoncini, Chaos 16, 023104 (2006). 



Periodic orbit bifurcations as an ionization mechanism 



20 




Figure 18. Relative ionization probability based on the maximum field rule, on 
Eq. (j4|), for (a) Case (I) with C = —2.42 and D — 2. Insets show the bichromatic field 
for = (bottom panel) and <fi — ir/3 (top panel). (b) Case (//) with C — —2.25 and 
D = 2. Insets show the bichromatic field for <f> — tt/6 (bottom panel) and 4> — 7r /2 
(top panel), (c) Case (III) with C = —0.24 and D = 0.3. Insets show the bichromatic 
field for (j> — ir/4 (bottom panel) and <f> — 37r/4 (top panel), and (d) Case (/V) with 
C = —0.625 and D = 0.55. Insets show the bichromatic field for cf> = (bottom panel) 
and <j) = 7r/5 (top panel). Arrows indicate the <p values for which the bichromatic fields 

arf drawn 



Periodic orbit bifurcations as an ionization mechanism 



21 



[25] J.G. Leopold and D. Richards, J. Phys. B 18, 3369 (1985). 

[26] P. Cvitanovic, R. Artuso, R. Mainieri, G. Tanner and G. Vattay, Chaos: Classical and Quantum, 

ChaosBook.org (Niels Bohr Institute, Copenhagen 2005). 
[27] J.M. Greene, J. Math. Phys. 20, 1173 (1979). 
[28] R.S. MacKay, Nonlinearity 5, 161 (1992). 
[29] J.R. Cary and J.D. Hanson, Phys. Fluids 29, 2464 (1986). 
[30] K.J. Schafcr and K.C. Kulander, Phys. Rev. A 45, 8026 (1992). 

[31] B.E. Sauer, S. Yoakum, L. Moorman, P.M. Koch, D. Richards, and PA. Dando, Phys. Rev. Lett. 
68, 468 (1992). 



